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diffusion along its field lines. We briefiy review the observational foundations of the Jansson- 
Farrar 2012 model for the large scale structure of the GMF, underscoring the robust evidence for 
a N-to-S directed, spiraling halo field. New results on the tensing effect of the GMF on UHECRs 
are presented, displaying multiple images and dramatic magnification and demagnification that 
varies with source direction and CR rigidity. 
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1. The Galactic Magnetic Field 

Our understanding of the Galactic magnetic field (GMF) has improved tremendously in recent 
years. The Jansson-Farrar (2012) (JF12) GMF model is currently the most realistic and compre¬ 
hensive model available [1,2]. In addition to the coherent component, the JF12 model describes 
the spatial variation of the random field strength and additionally has a striated random component, 
in both disk and halo. “Striated” is the term coined by JF12 to describe a field component that is 
locally aligned along some direction but with a randomly changing sign, as would originate from a 
random field being stretched or compressed in some preferred direction, e.g., by a wind emanating 
from the Galactic plane. 

The coherent and striated components of the JF12 model were constrained by fitting Faraday 
Rotation Measures of 40k all-sky extragalactic sources, simultaneously with WMAP polarized 
(Q,U) synchrotron emission maps [1]. They are proportional to the line-of-sight integrals of 
for the RMs and of for Q & U, where Hq and /icre are the thermal and cosmic ray electron 

densities. Fitting the WMAP total synchrotron emission map, allows the random field to be con¬ 
strained as well [2]. Altogether, there are more than 10,000 independent datapoints consisting of 
the mean RM, Q, U and I in 13.4 degree-squared Healpix pixels, each with its own astrophysical 
variance measured from the variance among the 16 sub-pixel measurements. The variance comes 
from Galactic random fields, local inhomogeneities in ^e,cre^ and in the case of RMs from the intrin¬ 
sic RMs and contributions of the extragalactic path. The robustness of the JF12 fit to using different 
masks and subtracting or not subtracting Galactic foreground where known, is a testimonial to the 
basic approach of de-weighting datapoints with large variance. See [1] for details. 

JF12 found their best fit when the orientation direction of the striated field is locally aligned 
with the direction of the coherent field, so that the striated contribution amounts simply to an 
enhancement of the polarized synchrotron signal relative to the RM signal [1, 2]. This suggests 
that compression of a coherent field, e.g., by SNe explosions, may be the origin of the striated field 
effect [3]. 

A major source of uncertainty in the GMF comes from our imperfect knowledge of the and 
ficre distributions, on which the observables depend. The ne,cre densities used by JF12 were from 
NE2001 [4] and GALPROP (A. Strong, private communication), but the GALPROP model is 
azimuthally symmetric, whereas a more realistic model is likely to have more structure, e.g., [5]. 
Fortunately, modeling the Galactic CR distribution is a topic of intense effort and more realistic 
models are beginning to emerge. The qualitative impact of a more structured fZcre is to reduce the 
strength of the JF12 field in the plane, as discussed in [3]. 

In spite of inadequacies in current models of AZg^cre^ the observational data on RM, Q, U unam¬ 
biguously demand a large scale opening-spiral component to the halo field with field lines running 
from South to North and the toroidal component having opposite signs in the N and S hemispheres. 
The relative signs of the poloidal and two toroidal components are consistent with what would be 
produced by differential shear in the Galaxy, starting with the poloidal component. Fig. 1 shows 
the field lines of the coherent halo field in the JF12 model; the field of the disk is rendered with finer 
but more dense lines in a different color and without directional arrows, to aid visual clarity. Fig. 
2 shows the RM, Q and U data (top line) and the predictions from the JF12 coherent+striated field 
model (2nd line; n.b.,the random field does not contribute to these observables). For comparison. 
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the predictions of the Sun and Reich (SRIO) [6] and Pshirkov, Tinyakov and Kronberg (PTKll) 

[7] models are also shown. The characteristic large scale L-R, N-S symmetry and antisymmetry 
pattern seen in the Q,U data calls for an outwardly-spiraling halo field as is present in the JF12 
model but not in the others. The Q,U data alone could be fit with a purely striated halo field, but 
the RM data demands that the halo-field has a coherent component of comparable magnitude to the 
striated component and determines its direction. Since and ficre are positive definite, modifying 
them affects the deduced strength rather than sign of the fields and cannot produce the pattern of 
the signs of RM, Q, U (barring an unnatural conspiracy between contributions from distant portions 
of the line of sight where the sign of the field can change). 

2. CR deflections in the GMF 

Using the NASA supercomputer Pleiades, and the Runge-Kutta CR propagation code CRT 

[8] , we have performed all-sky, fine-grained (1.26 x 10^ pixels) UHECR backracking through the 



Figure 1: The outwardly-spiraling halo field of the Galaxy (JF12); the field lines run from S to N. Visual¬ 
ization produced by T. Sandstrom, NASA. 
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Figure 2: Skymaps in Mollweide projection of (left to right) RM (in rad/m^), Stokes Q and Stokes U 
(in mK); Galactic longitude is / = 0° in the center and increases to the left. From top to bottom: row 1, 
data (white pixels (RM) lack data.); rows 2-4 simulated data from the JF12, SRIO and PTKll models, 
respectively. The 2, U data were not fit in a masked region where foreground contamination is greatest, see 
yet the JF12 fit is quite good in all regions; see [1] for the masks used. 

GMF, for rigitides of V, with IgR = 18.0, 18.1, ... 19.0, 19.2, 19.4, 19.5, 19.6,19.8 and 20. 
Inverting these maps, gives us forward-tracking from any specified extragalactic direction. We 
have also performed coarser-grained simulations (typically 10-40k pixels) for lower rigidities and 
recording the CR position every 30 yr, for a local source and a source at the Galactic Center, at 
rigidities IgR = 16.5, 17, 17.5, 18 and 18.5 to study the anisotropy of transient Galactic sources. 
Except otherwise noted, the plots below show results for the JF12 coherent field plus a specific 
realization of the random field obtained by weighting a Kolmogorov random field (KRF) having 
coherence length 100 pc, with the JF12 model for the spatial dependence of ^^and Some studies 
with other realizations of the KRF and with other coherence lengths have been performed [9, 10] 
and comparative plots are shown below as space permits. 

2.1 UHECR deflections 

Fig. 3 shows the observed arrival directions for 12 selected source directions - the centers 
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Figure 3: Arrival directions of 10 EV CRs, for 12 representative source directions marked with crosses, at b 
= ±41.5°, 1 = ±45° and ±135°, and at b=0,1 = 0, ±90° and 180°; events from a given source have the same 
color as the source marker. 


of the 12 base-level Healpix pixels - at 10 EV. One sees that deflections are generally very large. 
Except for the 4 directions farthest from the GC (b = ±41.5°, 1 = ±135°, the arrival directions 
are typically highly dispersed. The Galactic plane is something of an “attractor”, apart from the 
directions in the south and away from the galactic center, for which the deflections cause the arrival 
directions to come from further to the S. By comparing different realizations of the random field, 
one sees - not surprisingly - that the average deflection is basically the same as it depends mainly 
on the coherent field, while the degree of dispersion decreases when the coherence length or 5rand 
is reduced. Space does not permit inclusion of more examples, but as would be expected, mean 
deflections generally increase with decreasing rigidity. At higher rigidity, such that the deflections 
do not carry the UHECR too close to the Galactic plane, the deflections can be much smaller, 
especially for directions away from the Galactic center. 

The GME acts as a lens for UHECRs [11]. Although the total number of UHECRs is con¬ 
served, and an isotropic distribution of sources leads to an isotropic observed sky, the total flux 
from any given source direction can be magnified or demagnified relative to the flux in the absence 
of the GME, because smaller or larger areas of the plane wave of UHECRs coming from a given 
source direction can be focused onto the Earth (for illustrative figures see [3]). In the case that 
multiple distinct regions of the source plane are focussed onto Earth, multiple images of the source 
are produced, as can be seen in Eig. 3 and discussed in [9] for Cen A as the source. 

While the existence of magnetic lensing is to be expected, the actual impact of the Milky Way 
magnetic field on UHECR arrival directions can only now be explored, thanks to the new high- 
resolution simulations and use of a realistic field model. The results are quite remarkable. Eig. 4 
shows skymaps of logio of the magnification at various rigidities, as a function of source direction 
on the sky. Dark blue represents complete blindness to the given direction - UHECRs from those 
directions are deflected away by the GME and totally miss Earth - while the darkest red represents 
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Figure 4: Skymaps showing logic of tho magnification, for 4 different rigidities. Blue means no CRs from 
that source direction reach Earth (they are defiected away from us by the GMF) while a source direction 
indicated with the deepest red has a factor 30 magnification. For some source directions there is a rapid shift 
from magnification to demagnification as the rigidity changes. 

a magnification by a factor 30. As the random field is changed, e.g., reduced in strength or the 
coherence length changed, the details of the pattern of magnification changes, but qualitatively the 
picture is the same. The most important features of the magnification map are: 

• At high rigidity, the angular size of high-magnification and blind regions is small and their 
position varies with random field model. 

• For a given random field configuration, the magnification and demagnification for a given 
source direction can vary rapidly with rigidity. 

• As the rigidity decreases, the blind region grows until about half the sky is invisible - from 
behind the Galactic center and especially from the South. 

Fig. 5 (left) shows the trajectories of 20 randomly selected events of rigidity = E/Ze = 3 EV 
arriving to Earth through the coherent field only, or (right) with random field also, shortly before 
their arrival. 

2.2 VHE Galactic CR explosions 

Eig. 6 shows a snapshot of 0.3 EV CRs from an explosion at the Galactic Center, 100 kyr after 
the explosion. The important point to notice is how anisotropic the “diffusion” is - preferentially 
escaping vertically. Almost no CRs reach the solar circle at these energies. These results under¬ 
line the crucial importance of including anisotropic diffusion in a realistic magnetic field, before 
drawing conclusions about galactic cosmic rays, at any energy. 
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Figure 5: The tracks of 20 UHECRs chosen at random from an isotropically-arriving distribution, for rigid¬ 
ity = E/Ze = 3 EV; the left panel is with the coherent field only, and the right with a 100 pc random field. The 
magnification of sources away from and N of the GC, and blindness to sources behind the GC and especially 
to the south, is clearly seen. 
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Figure 6: The tracks of high energy (rigidity = E/Z = 0.3 EV) Galactic Cosmic Rays produced at the Galactic 
Center. Note the strong anisotropy in diffusion - almost no cosmic rays reach Earth before escaping the 
Galaxy. 
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